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^ by angular momentum loss (AML) driven by magnetic braking and gravitational 
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ABSTRACT 



radiation (GR) above the period gap, and solely by GR below the period gap. The 



O 

in, . . 

. best-fit revised model of CV evolution recently by lKnigge et al.l ( 120111 ). however 

indicates that AML rate below the period gap is 2.47(±0.22) times the GR rate, 
suggesting the existence of some other AML mechanisms. We consider several 
^ ! kinds of consequential AML mechanisms often invoked in the literature: isotropic 

lO wind from the accreting white dwarfs, outflows from the Langrangian points, and 

^ ! the formation of a circumbinary disk. We found that neither isotropic wind from 

the white dwarf nor outflow from the Li point can explain the extra AML rate, 
while onflow from the L2 point or a circumbinary disk can effectively extract the 
angular momentum provided that ~ (15 — 45)% of the transferred mass is lost 
from the binary. A more promising mechanism is a circumbinary disk exerting 
^ , gravitational torque on the binary. In this case the mass loss fraction can be as 

H : low as < 10^1 

Subject headings: stars: cataclysmic variables — stars: evolution — CVs: general 



1. Introduction 



Cataclysmic variables (CVs) are short-period binaries consisting of a white dwarf star 
accreting mat erial from a lower-mass main-sequence star that is overflowing its Roche-lobe 
( lWarnerlll995l ). In the standard model, mass transfer in CVs is driven by angular momentum 
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loss (A ML) driven by gravitatio nal radiation (GR) (IKraft et al.lll962l ) and magnetic braking 
(MB) JVerbunt fc Zwaanlll98lh . 



The obital period (Porb) distribution of CVs has been summarized by iRitter fc Klob 

(120031 ). which is bimodal with ~ 45% of CVs having period in the range of ~ 3 — 16 h, 
another ~ 45% with ~ 80 min —2 h, and the rest ~ 10% with ~ 2 — 3 h. The number of 
CVs in the period interval ~ 2 — 3 h is small, and it is known as the perio d gap. In seeking 



a pla u sible explanation of the period gap in CV evo lution, several authors (IRappaport et al. 



19831 : ISpruit k Ritterl Il983l : iLivio fc Pringld 119941 ) proposed the "disrupted" MB models 
which are related to the transition of AML mechanisms. The general ideal is that, mass 
transfer in CVs above the period gap is primarily driven by MB at a rate rapid enough, 
making the secondary star out of thermal equilibrium, so that the secondary star is oversized 
and has larger radius than a main-sequence star with the same mass. Along with mass 
transfer, the secondary star loses mass gradually, and finally become fully convective when 
Poi-b 3 h. AML is now assumed to been caused solely by GR, because MB is vanished. 
Accordingly the oversized secondary star begins to shrink and underfill its Roche-lobe in 
attempt to reach thermal equilibrium, cutting off matter transfer. CVs becomes very faint, 
and virtually unobservable. When the Roche-lobe is filled again due to orbit shrinking driven 
by GR at Porb — 2 h, mass transfer restarts from the secondary star. 



Most recently iKnigge et al.l (120111 ) reconstructed the complete evolutionary path fol- 
lowed by CVs, using the observed mass-radius relationship of their seco ndary stars. Fo r 



AML, they adopted a scaled versions of the standard GR loss rate and the IRappaport et al. 



(119831 ) MB law. With suitable normalization parameters, fan and /m_b, these recipes pro- 
vide acceptable matches to the observed data. The best-fitting scaling factors are fcR = 
2.47 ± 0.22 below the gap, and /mb = 0.66 ± 0.05 above, which describe the mass-radius 
data signicantly better than the standard model (/gr = /a/b = !)• 

Here we focus on the or igin of the enha.nced AML beloAV the gap, which has a l ready 



been mentioned befo re (e.g. iKolb fc Baraffd Il999l : iPattersonI l200ll : ISpruit &: TaamI 12001 



Barker fc Kolbl l2003l ) . One obvious candidate is the residual MB. Generally, the magnetic 
field of a low-mass star is supposed to be anchored at the interface between the convective 
envelope and the radiative core. For a fully convective star, such as the secondary star of CVs 
below the period gap, the interface disappears, and MB is thought to be closed. However, 
there is strong evidence that fully convective stars a re capable of generating significant 
magnetic fields (jSills et al.ll2000l : lAndronov et al.ll2003l ). which mig ht develop in a diffe rent 
way than low-mass main-sequence stars (for a discussion see §8.5 in iKnigge et al.ll201ll . and 
references therein). At present it is not clear whether the magnetically channelled stellar 
winds can produce MB strong enough to be consistent with the model for CV evolution (e.g. 
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Li et al.lll994[ ). In this paper we will examine other mechanisms that are needed to account 
for the extra AML. 

During the mass transfer processes, part of the transferred mass from the secondary 
star may escape from the binary system, carrying away the orbital angular momentum. 
AML associated wit h mass transfer is called consequential angular momentum loss (CAML) 
(iKing &: Kolbl 119951). There are several t ypes of CAML due to different way of mass loss in 



the CV evolution ( ISoberman et a 



and surroundin g accretion disk ( 



19971): fl) isotrop ic wind from the accreting white dwarf 



ang fc Kolblll995h . (2) outflow through the Lagrangian 



points Li or L9 (jVanbeveren et al 



I994J : iTaam fc Spruitll200lh . 



19981 ). and (3) a circumbinary (CB) disk ( van den Heuvel 



The organization of this paper is as follows. We introduce the basics of orbital evolution 
and possible CAML mechanisms in Section 2. In Section 3, we constrain possible AML 
machanisms to explain the extra AML rate below the period gap. Discussions and conclusions 
are presented in the final Section. 



2. Orbital evolution and CAML in CVs 

We first derive the r elatio n betw een the CAML ra t e and the mass transfer rate M2, 
following the approach in iKingI (119881 ) and iKnigge et al.l (120111 ) . 



We assume that the total angular momentum of the binary is dominated by the orbital 
angular momentum. 



J = M,M2{— 

where a is the binary separation. Mi, M2, and M are the masses of the white dwarf, the sec- 
ondary star, and the binary system, respectively. The relation between the binary separation 
a and the orbital angular velocity u is given by the Kepler's third law. 



GM 

Logarithmic differentiation of Eq. (1) gives 



J 



M 



_ Ml M2 a 
J~Mi^M2^2^~ 2M' 



(2) 



(3) 



The obital AML rate in CVs contains two items. One is the systemic AML due to MB 
and/or GR; the other is CAML, which is related to the mass transfer rate. So the total 
obital AML rate can be writen as 



J — Jsys + JcAML- 



(4) 
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When mass is transferred from the donor star to the white dwarf, we assume that a fraction 
6 of the matter flow escapes from the CV binary, i.e., 



M = 6M2, 

which means that the actual accretion rate of the white dwarf is 



(5) 



Ml = (5 - 1)M2. (6) 
Considering that the CAML rate is related to the mass transfer rate M2, we can write 



JcAML 



(7) 



where 1/ is a para meter as a function of 5 ( iKnigge et al.l 1201 ll ) . In order to eliminate d/a in 
Eq. (3), we use the iPaczynskil (jl97ll ) formula for the Roche-lobe ridius Rl of the donor star. 



Rl = 0A62{M2/My/^a. 



Logarithmic differentiation of Eq. (8) yields 



Rl 
Rl 



M2 M a 
3M2 ~ 3M^a' 



(8) 



(9) 



For a Roche- lobe-filling donor star with steady mass transfer, variation of the stellar radius 
i?2 and the Roche- lobe radius Rl should be in step, i.e.. 



Rl 
Rl 



R2 
R2 



(10) 



To deal with i?2, we use the mass-radius relationi?2 = M2 and its logarithmic differentiation 



R2 
R2 



c 



_M2 

Mo' 



(11) 



where ( is the mass-radius exponent. iKnigge et al.l (I2OIII ) showed that after a CV donor 
emerges from the bottom of the period gap, the exponent ( evolves from 0.8 (in thermal 
equil ibrium) to ~ 1/3 at the minimum period, and finally to ~ —1/3 (see also lRappaport et al, 



19821 ). They adopted broken-power-law fit to the updated M2 — R2 data in iKniggd (120061 ). 
with C = 0.3 for M2 < 0.07 M©, and C = 0.61 for O.OTMq < M2 < 0.2 Mq. Here O.2M0 
and O.O7M0 are the secondary masses just below the gap and at the minimum period where 
the secondary stars become degenerate and the orbital periods start to bounce back into 
the period- increasing phase, respectively. Adopt q = M2/M1 and typical white dwarf mass 
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Ml = O.QMq in CVs, C decreases from ~ 0.61 to ~ 0.3 with decreasing q. In the following 
we calculate ( by using the M{R) relation derived fro m numerical calcula tions of CV binary 
evolution, which is very close to the empirical one in lKnigge et al.l (l201ll ). 



Combining Eqs. (3)- (11), one can derive 

M2 _ Jsys 

W2 ~ Id' 



where 



Mo Mo 



3M' 



u. 



(12) 
(13) 



Here the systemic AML mechanism below the gap is GR, and its rate is given by (iLandau &: Lifshitz 

32 



19511 ) 



Jgr 
J 



5 c5 



MiMaa" 



(14) 



Next we consider the possible CAML mechanisms. For isotropic wind (e.g., outflows 
due to nova explosions and/or wind emanating from the accretion disk), we assume that it 
carries the white dwarf's specific obital angular mometum ji, 

M2 



Ji 



:J. 



MiM 

For outflows from the Lagrangian point L2 the specific obital angular mometum is 



32 



(15) 



(16) 



where a^^ is the distance between the mass center of binary and the L2 point. Finally for 
CB disks we assume that it extracts the orbital angular momentum from the binary by tidal 
force. The specific obital angular momentum is given by. 



J3 



(17) 



wher e 7^0 is taken as the inner radius of the CB disk. Usually, 7 = 1.5 (ISoberman et al. 
1997h . 



3. Constraining the CAML machanisms below the period gap 



The best-fit revised model of CV evolution in iKnigge et al.l (120 111 ) indicates that the 
AML rate below the period gap is 2. 47 Jgr- This means that some other AML machanisms 
besides GR should work. In the following we will discuss the feasibility of the AML related 
to isotropic wind, outflow from L2 point, and CB disks, respectively. 
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3.1. Isotropic wind 

Combining Eqs. (5), (7) and (15), we can derive tlie obital AML rate, Jcaml,i of 
isotropic wind as 

and hence 

MiM ^ ^ 

Let JcAML,i = l-47JGij and combined Eqs. (12)-(14) and(18)-(20), we have 

2.47 5M| , , , , Ma ,,M2 Ms, 

which can be transformed into 

(2,5 + 1.5C-39)(l + g) 

For CVs below the period gap, < g < 0.33, C > 0, so 5 is always negative, indicating that 
the isotropic wind cannot provide the extra 1A7Jgr AM L. This conclusion re mains valid if 



we let most of the material to escape from the Li point (IBarker fc Kolbl 120031 ). 



3.2. Outflow from the L2 point 

Similar as in section 3.1, we derive the obital AML rate, Jcaml,2, of outflow from the 
L2 point by combining Eqs. (5) and (16), 

jcAML,2 ^ A^M2 .22) 

J a2 Ml M2 ' ^ ^ 

which gives 



Then we obtain the following relation 



L47^^ = + - Ml + ^^M: - 3m)- 
After simplification we have 

^ (2.5 + 1.5C-3g)(l + g) 

" 5(l + g)2x2-g(2 + 3g)' ^ ' 



where x = ai^/a. For different g, the values of x are given in lMochnackil ( 119841 ). The relation 
between 5 and q is shown in Figure 1. It is seen that the required value of 5 ranges from 
~ 0.15 to ~ 0.45. 
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3.3. CB disks 



The origin of tlie CB disks may stem from the remnant of the late stage of the common 
envelope evolution phase that formed the CV s, or matter outflow from CVs during the mass 
transfer processes (ITaam fc Squdquistll2000l ). Through tidal interaction of the CB disk and 
CVs binary, the CB disk can extract obital angular momentun from binary if part of the 
transferred mass flows into the disk rather onto the white dwarf. The corresponding AML 
rate is, 

jcAML,3 = S-fM2a^u. (26) 

Then we can derive 

. (2.5 + 1.5C-3g)(l + g) 



7.56(1 + g)2 -g(2 + 3g)' 
The relation between 6 and q is shown in Figure 2. The value of 6 lies between 
~ 0.4, comparable with that in the case of outflow from the L2 point. 



(27) 
0.22 and 



The above result is under the assumption that AML is only caused by mass loss. Ac- 
tually there is gravitational torque between the CB disk after it is f ormed and the binary , 
which is more efficient to extract angular momentum from the binary (jSpruit &: Taanj|200ll ). 
The AML rate under this torque can be expressed as 



J, 



CB 



^(__)5M2( 



,1/3 



orb 



{21 



where t is the time since mass transfer begins, and t^i is the viscous timescale at the inner 
edge of the CB disk, give n by = 27^Porb/3vra;/3^, where a is the viscosity parameter 
(IShakura fc SunyaevI 119731 ) . and (3 is the ratio of the scale height to the radius of the disk. 
Equation (28) can be further simplified to be 

JcB = A{GMf/^6M2t^^^ (29) 

where A = (3a/3^/4)^/^. Similar as the derivation of Eq. (20), we obtain the following relation 
for the CB disk, 

2.47 AGi/e^vf 7/6^^1/3 ^^^^ 



1.47 MiaV2 



Mo Mo 
(5/6 + (/2)-^ + .(^ 



3M^ 



(30) 



Previous investigations (ITaam fc Spruitll200ll : iTaam et al.ll2003l : IWillems et al.ll2005l ) suggest 
that very small values of (5(<C 1) are required for CV evolution. Thus the third term on the 
rhs of Eq. (30) can be neglected. Considering the fact that C/2 and q also roughly counteract 
each other, and M ~ Mi, Eq. (30) is changed to be 



orb X 1/3 



2A(27r)i/3^ t 
X lO^^an of /3( 



*0.01 1^0.03 



-2/3^- Porb,2 xi/3 



(31) 
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where ao.oi = a/0.01, /3o.o3 = 
This is close to the result 6 



/3/0.03 flBelle et al.ll2004f ). Porb,2 = Porb/2 hr, and t g = t/10^ yr 
3 X 10~^ adopted in the numerical calculations by iTaam et al. 



([20031) for CV evolution. 



To investigate the effect on AML of the uncertainties in treatment of CB disks we have 
numerically solved Eq. (30). Figures 3 and 4 show the calculated values of 5 as a function of 
q for different values of a and /3. It is seen that generally smaller a or P corresponds to larger 
5. This is easy to understand with Eq. (29): smaller a or /3 indicates less efficient angular 
momentum transfer within the disk, which requires more mass input into the CB disk to 
guarantee enough AML from the orbit. Nevertheless, the values of 6 are always small (< a 
few 10~^) when we change a from 0.001 to 0.1, and (3 from 0.005 to 0.1. This implies that CB 
disks are indeed very efficient in draining orbital angular momentum through gravitational 
torques even with a very small mass input rate. 

To show how CB disks can influence the evolution of CVs below the gap, we have also 
performed binary evolution c a lculations a dopting ar i updated version of the stellar evol ution 
code developed bv lEgdetonl Jl97ll . Il972h (see also lflan et al~lll994l : Fols et al. Ill995[ l. We 
set initial solor chemical compositions (i.e., X = 0.7, Y = 0.28, and Z = 0.02) for the donor 
star, and take the ratio of mixing length to the pressure scale height to be 2.0, and the 
convective overshooting parameter to be 0.12. 

We follow the evolution of a CV just below the gap with a donor star of mass 0.2Mq 
and an orbital period ~ 0.1 d. For CB disks we take A6 as one free parameter to assess its 
influence on the evolution of CVs, since a, /3 and 6 are always combined together i n evaluating 
the A ML rate (see Eq. [29]). For typical values of a{= 0.01) and (3{= 0.03) ( iTaam et al. 
20031 . and references therein), A ^ 0.02. Considering the fact that suita ble value of S may 
range from ~ 10~^ to a few 10"'* ( iTaam et al.l 120031 : IWillems et al.l 120051 ). we constrain the 
adopted value of AS to be less than ~ 10~^. Larger AS may cause unstable mass transfer. 
The stability of mass trans fer can be examiiied by comparing the Roche-lobe radius exponent 
(l due to mass loss with ( (jSoberman et al.lll997l ). which are shown in Figure 5 as a function 



of q. From top to bottom, the solid curves represent Cl with AS ranging from 5 x 10^^ to 
1 X 10~^ in steps of 1 x 10~^, and the dashed lines show the mass-radius exponent ( of 0.8, 
1/3, and —1/3, respectively. Since the mass transfer would be unstable when > C; we can 
derive that AS should be < (2 — 3) x 10~^ to guarantee stable mass transfer. 

In Figure 6 we present examples of the evolutionary sequences of the donor mass, orbital 
period, mass transfer rate, and the ratio of total AML rate and the AML rate due to GR 
(A^ = J/ Jgr) to show the influence CB disks. The panels from top to bottom correspond to 
AS = (i.e., AML due to GR sololy), 4 x 10~^ 9 x 10"^ and 2 x 10-^ respectively. When 
AS = Ax 10~^ (or 5 = 2 X 10"^(y4/0.02)~^), the evolution seems similar to that with GR only. 
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However, the mass transfer processes are actually accelerated, with ~ 10% higher average 
mass transfer rate and ~ (10 — 15)% higher AML rate than in the GR-only case. This 
tendency becomes more intense when A6 goes up. When A6 = 2 x 10~^, the average mass 
transfer rate is enhanced by ~ 70%, and the AML rate J becomes ~ 2 — 3 times Jgr- These 
results imply that a tiny fractional input rate [S < 10"'^) into CB disks can significantly 
change the evolution. This is in contrast with mass loss-assissted AML mechanisms, which 
usually require a much larger fraction of the transferred mass to leave the binary system 
(although in real situation there might be multiple mechanisms work simultaneously). 



4. Discussion and Conclusions 



The secular evolution of CVs is thought to be driven by the AML. Two mechanisms 
usually invoked to account for the dissipation of AM are GR and MB of the secondary 
star. However, this dual loss-mechanism cannot completely account for the magnitudes of 
the mass-transfer rates inferred for some CV s, and for the large sp read in the mass transfer 



rates observed at a given orbital period (e.g. ISpruit fc Taamll200ll ). Likely additional AML 



mechanisms include mass loss from the binary during the mass transfer process, which carries 
away the orbital angular momentum from the binary. Below the period gap in CVs evolution, 
AML mechanism was usually considered to be driven by sol ely GR, while the b est-fit result 
with observations indicates that AML rate is about 2.47 Jgr (iKnigge et al.ll201ll ). This offers 
a possibility to constrain the AML mechanisms besides GR. We consider several kinds of 
CAML mechanisms often invoked in the literature: isotropic wind from the accreting white 
dwarfs, outflows from the Langrangian points, and the formation of a CB disk. 

We found that neither isotropic wind from the white dwarf nor outflow from the Li 
point can explain the extra 1A7Jgr AML rate, while outflow from the L2 point or a CB disk 
is more effective in extracting the angular momentum. For a O.GMq white dwarf, the fraction 
6 of mass loss in the total transferred mass is ~ 0.15 — 0.45 or ~ 0.2 — 0.40, respectively. 
Actually it is found that 6 is allways less than 0.45 for different mass of the white dwarf in 
our calculations. Note that when mass is lost fro m the L2 point, it is very likely to form a 
CB disk around the binary (Ivan den Heuvellll994l ). so it is not surprised that the values of 6 
are very close in these two cases. 

When the tidal interaction between the CB disk and the binary is included, the mass 
transfer rate can be enhanced much more efficiently, and a very small fraction {6 < 10~^) 
of mass loss is required. This will suggest a much lighter CB disk than in the former cases. 
The CB disks are thought to be large (up to several AU in radius) and cool (a few thousand 
K at the inner edge to less than 1000 K at the outer edge), with peak emission in the 
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infrared ( iDubus et al.ll2002l ). Although detection of excess infrared emission from magnetic 
CVs provides obs ervational support for the presence of cool gas (and possibly a CB disk) 



surro unding CVs (IHowell et al.ll2006l : iDubus et al.ll2007l : iBrinkworth et al.ll2007l : iHoard et al. 



20071 ). there are still many open questions associated with the formation of the CB disks. 
Future observations with the measurement of the disk masses might distinguish the ways of 
angular momentum transfer between the CB disk and the CV binary. 
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by the Natural Science Foundation of China (under grant numbers 10873008 and 11133001), 
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Fig. 1. — The relation between 6 and q in the case of outflow though the L2 point. 
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Fig. 2. — The relation between S and q in the case of CB disks. 
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Fig. 3. — The fractional mass input rate, 5, as a function of q for CB disks with gravitational 
interaction with the binary. Here we take /3 = 0.03, and the curves from top to bottom 
correspond to a; = 0.001, 0.005, 0.01, 0.05, and 0.1, respectively. 
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Fig. 4. — The fractional mass input rate, 5, as a function of q for CB disks with gravitational 
interaction with the binary. Here we take a = 0.01, and the curves from top to bottom 
correspond to /3 = 0.005, 0.01, 0.03, 0.05, and 0.1, respectively. 
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Fig. 5. — The Roche-lobe radius exponent as a function of q. The five sohd curves represent 
C,L with different A5 from 5 x 10~^ to 1 x 10~^ in steps of 1 x 10~^ (from top to bottom). 
The dashed hues from top to bottom correspond to the mass-radius exponent Q = 0.8, 1/3, 
and -1/3, respectively. 
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Fig. 6. — Evolution of the donor mass (Mq), orbital period (d), mass transfer rate {MqJt'^), 
and the ratio (A*") of total AML rate and the AML rate due to GR. The panels from top to 
bottom correspond to AML due to GR sololy, and both GR and CB disks with AS — Ax 10~^, 
9 X 10"®, and 2 x 10"^, respectively. 



